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i The perturbative renormalization group(RG) equation is applied to resum divergent 

. series of perturbative wave functions of quantum anharmonic oscillator. It is found that 

the resummed series gives the cumulant of the naive perturbation series. It is shown that 
a reorganization of the resummed series reproduce the correct asymptotic form of the wave 
function at x — * oo when the perturbation expansion is stopped at the fourth order. A 
brief comment is given on the relation between the present method and the delta-expansion 
method, which is based on a kind of a nonperturbative RG equation. 



§1. Introduction 



> 

^0 ■ In this talk, we shall present a novel use of the perturbative renormalization 

as : group equations a la Gell-Mann-LowS) to solve problems in other fields than quantum 

field theory and statistical mechanics. 

When applying a perturbation theory to any problem, one will encounter diver- 
gent or at best asymptotic series. One needs to resum the divergent series to obtain 
a sensible result from the perturbation theory. Thus various resummation tech- 
niques have been devised.™ Recently, a unified and mechanical method for global 
and asymptotic analysis has been proposed by Goldenfeld et al.cP This is called 
the renormalization group (RG) method. Indeed the RG equations a la Gell-Mann- 
Low is known to have a peculiar power to improve the global nature of functions 
obtained in the perturbation theory in quantum field theory© A unique feature of 
Goldenfeld et al's method is to start with the naive perturbative expansion and allow 
secular terms to appear in contrast with all previous methods; a adding unperturbed 
solutions to the perturbed solutions so that the secular terms vanish at a "renormal- 
ization point" t = to and then applying the RG equation, one obtains a resummed 
perturbation series. 

Subsequently, the RG method has been formulated on the basis of the classical 
theory of envelopes; a>>E> it has been indicated© that the notion of envelopes is also 
useful for improving the effective potential in quantum field theory. |) til 

The purpose of the present work is to make an introductory account of the RG 
method as formulated in Ref.'s ||),[5|) and apply the method to Schrodinger equation 
of the quantum anharmonic oscillator (AHO) to resum the naive perturbation series 
by the RG equation; we shall also examine whether our method is useful to obtain 



*' E-mail address: kunihiro@rins.ryukoku.ac.jp 
*rj Shirkov has extracted the notion of functional self-similarity (FSS) as the essence of the 
RG.EP We remark that the notion of FSS is only applicable to autonomous equations, while that of 
envelopes is equally applicable to non-autonomous equations, too. 



typeset using ■pTpJ^.sty <ver.0.8> 



2 



T. Kunihiro 



the asymptotic form of the wave functions. 

AHO has been a theoretical laboratory for examining the validity of various 
approximation techniques.^ Bender and BettencourtllH) have recently shown that 
multiple-scale perturbation theory (MSPT) can be successfully applied to the quan- 
tum anharmonic oscillator; see also Ref. |TT| ). Using the fourth order perturbation 
series, they obtained an asymptotic behavior of the wave function ip{x) for large x 
which is in good agreement with the WKB result. One should remark here that a 
further resummation had to be adopted for the latter case, which is not intrinsic in 
MSPT and a similar method had been proposed by Ginsburg and MontrollEJ). 

We shall show that the resummation of the perturbation series of the wave 
functionsis performed in the RG method more mechanically and explicitly than in 
MSPT: lH Since our method is easy to perform to higher orders of the perturbative 
expansion, it is examined how the results of Bender and Bettencourt persist or are 
modified at the higher orders. 

In the next section, an introduction is given of the RG method for improving per- 
turbative solutions of differential equations. Then the method is applied to quantum 
mechanics in §3. The final section is devoted to a summary and remarks. 

§2. The RG method for global analysis; simple examples 

In this section, using simple examples, we shall show how the RG method B 
works for obtaining global and asymptotic behavior of solutions of differential equa- 
tions. One will see that the reasoning for various steps in the procedure and the 
underlying picture are quite different from the original ones given in Ref. ^): The 
present formulation emphasizes the role of initial conditions and the relevance to 
envelopes of perturbative local solutions. ffl'H 

2.1. A linear equation 

Let us first take the following simplest example, B 

d 2 x dx 

Hfi +e H + x = - {2A) 

where e is supposed to be small. The solution to (2.1) reads x(t) = ^4exp(— et/2) 
• sin(y / l — e 2 /4i + 9), where A and 9 are constants. 

Now, let us obtain the solution around the initial time t = to in a perturbative 
way, expanding 

x(t,t ) = x (Mo) + ezi(Mo) + e 2 x 2 (Mo) + (2.2) 
The initial condition may be specified by 

x(t ,t ) = W(t ). (2.3) 



*' As for the application of the RG method a la Goldenfeld et al, to the time-dependent 
Schrodinger equations, see Ref. p^[),p^). 
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We suppose that the initial value W(h) is always on an exact solution of Eq.(2.1) 
for any to. We also expand the initial value W(h); 

W(t ) = W (t ) + eWxih) + e 2 W 2 (h) + (2.4) 

and the terms Wi(h) will be determined so that the perturbative solutions around 
different initial times to have an envelope. Hence the initial value W(t) thus con- 
structed will give the (approximate but) global solution of the equation. 

Let us perform the above program. The lowest solution may be given by 

x (t, t ) = A(t ) sia(i + 9(h)), (2.5) 

where it has been made it explicit that the constants A and 9 may depend on the 
initial time to- The initial value W(h) as a function of to is specified as 

Wo(t ) = x (t , to) = A(t ) sin(t + 0(t )). (2.6) 

We remark that Eq.(2.5) is a neutrally stable solution; with the perturbation e ^ 
the constants A and 9 may move slowly. We shall see that the RG equation, which 
will be interpreted as the envelope equation, & gives the equations describing the 
slow motion of A and 9. 

The first order equation now reads x± + x\ = — Acos(t + 9), and we choose 
the solution in the following form, 

xi(t, to) = -A/2 ■ (t - to) sm(t + 9), (2.7) 

which means that the first order initial value W\(h) = so that the lowest order 
value Wo(*o) approximates the exact value as closely as possible. Similarly, the 
second order solution may be given by 

x 2 (t) = A/8 ■ {(t - t ) 2 sin(t + 9) - {t - t ) cos(t + 9)}, (2.8) 

thus 1^2(^0) = again for the present linear equation. 

It should be noted that the secular terms have appeared in the higher order 
terms, which are absent in the exact solution and invalidates the perturbation theory 
for t far from to- However, with the very existence of the secular terms, we could 
make Wi(to) (i = 1, 2) vanish and W{to) = Wo{to) up to the third order. 

Collecting the terms, we have 

x(t, t ) = A sin(i + 9) - eA/2 ■ (t - t ) sm(t + 9) 

+e 2 A/8 ■ {(t - t ) 2 sm(t + 9) - (t - t ) cos(i + 9)}, (2.9) 

and more importantly 

W(t Q ) = W (to) = A(t ) sin(t + 0(h)) (2-10) 

up to 0(e 3 ). Notice that W(h) describing the solution is parametrized by possibly 
slowly moving variable A(h) and 4>(to) = h + 9(h) in a definite way. 
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Now we have a family of curves {Ct }t given by functions {x(t, to)}t parametrized 
with to. They are all on the exact curve W(t) at t = to up to 0(e 3 ), but only valid 
locally for t near to- So it is conceivable that the envelope E of {Ct }t which con- 
tacts with each local solution at t = to will give a global solution. Thus the envelope 
function x E (t) coincides with W(t); 

x E (t) = x(t,t) = W(t). (2.11) 

Our task is actually to determine A(to) and 9 (to) as functions of to so that the family 
of the local solutions has an envelope. According to the classical theory of envelopes, 
the above program can be achieved by imposing that the envelope equation 

= 0, (2.12) 

dto 

gives the solution to = t.a Inserting (2.9) into (2.12), we have 

§ + eA = 0, £ + £=0, (2.13) 
dto dto » 

where we have utilized the fact that dA/dt is 0(e) and neglected the terms of 0(e 3 ). 
Solving the equations, we have 

A(t ) = Ae- et °/ 2 , 9(t ) = -e 2 /8-t + 9, (2.14) 

where A and 9 are constant numbers. Thus we obtain 

x E (t) =x(t,t) = W (t) = Aexp(-e/2 ■ t) sin((l - e 2 /8)t + 9), (2.15) 

up to 0(e 3 ). Noting that y/l - e 2 /4 = 1 - e 2 /8 + <3(e 4 ), one finds that the resultant 
envelope function x E (t) = Wo(t) is an approximate but global solution to Eq.(2.1). 

2.2. Forced Duffing Equation 

As an example of nonlinear equations, we take the forcedJOumng equation where 
the frequency of the external force is near the intrinsic one;@M 

x + 2e^yx + (1 + ea)x + ehx 3 = e/cost, (2.16a) 
y + 2e-yy + (1 + ea)y + ehy 3 = e/sint. (2.16b) 

When 7 = / = 0, the equation describe the usual (classical) anharmonic oscillator 
with the quartic potential. Defining a complex variable z = x + iy, one has 

z + 2e 7 i + (1 + ea)z + eh/2 • (3\z\ 2 z + z* 3 ) = efe u . (2.17) 

We suppose that e is small. 

Expanding z as z(t,to) = zo(t,to) + ez±(t,to) + e 2 Z2(t,to) + • ■ ■, we first, as 
before, obtain a perturbative solution around the initial time t = to, with the initial 
condition 

z(to,t ) = W(t ). (2.18) 



*' There are some errors in the Appendix of Ref. which are corrected below. 



Renormalization- group Resummation of the Perturbative Wave Functions 5 



The initial value W(io) is to be determined so as to the perturbative solutions at 
to varied continue smoothly. We thus also expand the initial value as W(to) = 
W (t ) + eW 1 (t ) + e 2 W 2 (t ) + ---. 

A simple manipulation gives the solution up in the first order approximation as 

z (t- 1 ) = A{t )e il - ie/2 • (t - t )[f - A{t ){a + 2i 7 ) 

-3/i/2 • \A{t )\ 2 A{to)]e u + e/32 • A*{t Q fe~ m + 0(e 2 ). (2.19) 

This solution implies that we have chosen the initial condition as 

W(t ) = z(t , t ) = A(t )e lt + e/32 • A*(t ) 3 e 3U + 0(e 2 ). (2.20) 

The RG or envelope equation reads 

dz 

— = 0, (2.21) 
ato 

with to = t, which leads to 

A = ie/2 ■ [(a + 2i-/)A + 3h/A ■ \A\ 2 A - /], (2.22) 

up to 0(e 2 ). Here we have discarded terms such as edA/dt, which is 0(e 2 ) because 
dA/dt = 0(e). The envelope is given 

z E (t) = z(t; t =t) = A(t)e lt + ^A* 3 e~ 3it + 0{e 2 ). (2.23) 

We identify ze(£) with a global solution of Eq.(2.17); x(t) = Ke[zE] and y(t) = Im[z^] 
are solutions to Eq.(2.16). If we parametrize the amplitude as A(t) = R(t)exp(i6(t)), 
we have 

x(t) = R(t) cos(t + 9 (t)) + e— R(t) 3 cos 3(t + 6 (t)), (2.24a) 

32 

y(t) = R(t)sm(t + e(t))-e^-R(t) 3 sm3(t + 6(t)), (2.24b) 
where the real functions R(t) and 0(T) satisfy the following coupled equation; 



R(t) = - ~ [y/a 2 + 47 2 sin aR(t) + fsin9 (t)] , (2.25a) 



R(t)9(t) = - [^a 2 + A~i 2 cos aR(t) + -hR(t) 6 - f cos 0(t)], (2.25b) 
with tana = 2j/o~. 

When the external force is absent, i.e., / = 0, Eq.(2.25) is readily solved, yielding 

R(t) = R e~^ 1 , 9(t) = Aujt - —RHe- 2 ^ - 1} + 6 , (2.26) 

87 



with 7 = 1/2 • \J g 2 + 47 2 sin a, Aoj = e\/o~ 2 + 47 2 cos a, Ro and #0 being constantM 
Inserting R(t) and 8(t) into x(t) and y(t) in (2.24), one has an explicit form of the 



**> When 7 = 0, R(t) = Ro = const, and 9(t) = e(a + 3hR 2 /4)t + o . 
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solution to the anharmonic oscillator. Notice that the solution contains an infinite 
order of e, although we have only made a first order calculation. When / ^ 0, the 
coupled equation can not be solved by quadrature hence the equation is not reduced 
to a simpler equation than the original ones, which may be a reflection that the 
forced Duffing equation gives chaotic solutions. 

Extensive applications of the RG method to nonlinear equations including partial 
differential equations and difference equations are given in Ref.'s p|),|4|),[T5|), 

§3. The RG-resummation of divergent perturbative wave functions 

Now let us apply the RG method to a quantum mechanical problem. We take 
the anharmonic oscillatorO 

(H - E)ip(x) = 0, H = p 2 + x 2 /A + ex 4 /4, (3.1) 

with the boundary condition ^(±oo) = 0. WKB analysis shows that for large x, 

i>(x) ~ exp{-^|x| 3 /6}. (3.2) 

We are interested in how the perturbation theory can reproduce the WKB behavior 
or not. 

3.1. Preliminaries 

We first apply the Bender- Wu method for performing Rayleigh-Schrodinger 
perturbative expansion: 



oo 



V>(x) ~ en Vn(?) and E(e) ~ ]T e n E n . (3.3) 

ra=0 ra=0 

The boundary condition is taken as yo(0) = 1 and y n >i(0) = 0. The lowest order 
solution reads yo(x) = e~ x 7/4 and Eq = 1/2. The higher order terms with n > 1 
are written as 

y n (x) = e- x2 / A P n (x), (3.4) 

where P n {x) is a polynomial, which satisfies the recursion relation; 

4 n— 1 

i*'(x) - xP^x) = ^-P n -i{x) - Pj(x)E n -j- (3.5) 

We remark that y n {x) in Eq.(3.4) may be identified as a secular term because it is a 
product of the unperturbed solution and a function that increases as x goes large. 
The eigenvalues E n are given by 

1 f°° x A n_1 

E n = —= dx y ( x ){—P n ^{x) - V P 3 (x)E n . 3 }. (3.6) 

Notice that E n is determined in terms of only the polynomials Pj(x) with j < n — 1. 

With these eigenvalues Ej (j < n), P n {x) is determined by Eq.(3.5). The 
polynomials P n (x) up to the six order are presented in Ref.10), which we refer to. 

*' We use the notations of Ref. 10) for later comparison. 
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3.2. The RG method 

Now we apply the renormalization group method as formulated in §2. 
We first calculate the wave function ip(x; xq) around an initial point x = xq in a 
perturbative way: 

oo oo 

ij(x;x ) ~ e>n(x;x ) and E(e) ~ ^ e n E n , (3.7) 

n=0 ra=0 

with the initial or boundary condition (BC) at x = xo; 

V>(x ;x ) = W(x ). (3.8) 

We suppose that the boundary value W(xq) is always on an exact solution of Eq.(3.1). 
We also assume that W{xq) can be also expanded in a power series of e 

oo 

W(x ) = t n W n (x ) (3.9) 

n=0 

If we stop at a finite order, say, the iV-th order, we will have J2n=o eU< Pn( x ) = 
(x; xq) which is valid only locally at x ~ xq. However, one may take the following 
geometrical point of view; namely, we have now a family of curves {ip^ (x;xq)} Xo 
parametrized with xq, and each curve of the family is a good approximation around 
x = xq. Then, if each curve is continued smoothly, the resultant curve will valid 
in a global domain of x. This is nothing else than to construct the envelope of 
the family of curves. More specifically, we only have to determine the boundary 
values W(xq) so that the perturbative solutions around x = xq form an envelope; 
see Fig.l. Furthermore, to be as accurate as possible, the lowest value Wq(xq) should 
approximate the exact value i(j(xo,xo) as close as possible, or W n >i(xo) should be 
made as small as possible. 



V( x ) W(x ) 




V|/(x,x' ) 



x x'o x 

Fig. 1. The local solutions and the global solution as their envelope. 
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The above program can be performed as follows. The decisive observation is 
that the lowest order solution may be written as 

Mx;x ) = Z(x )e- x2 /\ (3.10) 

with an arbitrary function Z(xq) of We remark that the Eq = 1/2 in this case, 
too. Accordingly, the boundary value is 

W (x ) = Z{x )e- x y\ (3.11) 

The higher order terms may be written as 

ip n (x;x ) = Z(x )e~ x2/A Q n (x;x ), (3.12) 

where Q n (x; xq) is a polynomial of x, dependent on xq. Q n (x; xq) is found to satisfy 
the same recursion relation Eq.(3.5) as P n (x). We impose the boundary condition 
(B.C.) as 

ip n (x ]x ) = W n (x ) = or Q n (x ;x ) = 0, (3.13) 

for n > 1 because we want to make the lowest order solution as close as the exact 
one at x = xqEIA Thus one readily obtains 

Q 1 (x;x ) = P 1 (x)-P 1 (x ),. (3.14) 

We remark that a constant is the solution of the homogeneous equation!***)! The 
second order equation reads 

4 2 

T 

-Pi(x )(x 4 /4-£i). (3.15) 



(Qfi''~ 
— Px{x) - 2J Pj(x)E n -j 

j=0 



One can verify that Ei is given by Eq.(3.6). Since Eq.(3.15) is a linear equation 
and the inhomogeneous part is a linear combination of those for Pi{x) and Pi(x), 
Q2{x;xq) satisfying (3.13) is given by 

Q 2 (x;x ) = (P 2 (x) - P 2 (x )) - P 1 (x )(P L (x) - Pi(x )). (3.16) 

Repeating the procedure, one finds that Q n (x;xo) are expressed in terms of Pj(x) 
(j < Ti); see Ref. ^) for the explicit forms for n < 6. 



*' It may be noted that the amplitude Z{xq) and the boundary position xq correspond to 
the wave function renormalization constant and the renormalization point in the field theory, 
respectively. 

**' We remark that the B.C. (3. 13) or Eq.(3.18) below corresponds to a renormalization condition 
in the field theory. 

***) w e remark that the second term of Eq.(3.14) corresponds to a counter term in the renormal- 
ization procedure in the field theory. 
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Thus we obtain the approximate solution valid around x ~ xo, 

oo 

V>(x; x ) ~ Z(x )e- x2 / 4 ^ e n Q n {x; x ), 

n=0 

where 

^(x ;x ) = W (x ) = Z(x Q )e- x V\ 



(3.17) 



(3.18) 



Now we have obtained a family of curves {Y>(x; xq)} xq with xq parametrizing the 
curves; each curve is a good approximation for x around x - The envelope of the 
local solutions is giving by solving the following equation, 



dtp(x;x ) 



dxo 



0, 



(3.19) 



XQ=X 



which is in the same form as the renormalization group equation; one may have 
called it the RG equation. The equation gives a condition which Z(xq) must satisfy; 



iry j OO 

u n=0 



Xq=X 



Defining f n (x) by 



one obtains 



--^—Qn(x;x ) 
dxo 



XQ=X 



dfnjx) 

dx ' 



Z(x) = A-exp[£ e n f n (x)], 



(3.20) 



(3.21) 



(3.22) 



n=0 



where A is a constant. With this solution, the global solution i/je(x) is given by the 
boundary value by construction; 



i, E {x) = W (x) = Z(x) e ^ 2 / 4 = A ■ exp[-^ + ]T e n f n (x)}. (3.23) 



This is one of the main results of the present work. 
/ n (x)'s are easily calculated in terms of P n (x); 



h{x) =Pi(x) = --x 2 -— - 4 



1 



n=l 



21 



, - — x 4 ,/ 2 (x) = P 2 (x) - ^Pi(x) 2 = ^x 2 + -^x 4 + ^x 6 , 

64 



256 ' 



2 v ' 16 

/ 3 (X) = P 3 (X) - P(x)P 2 (x) + 1P 1( X) 3 = -^X 2 - |x 4 - |L X 6 _ 

/ 4 (ar) = P 4 {x) - Pi(x)P 2 {x) - ^ 2 (*) 2 + P(x) 2 P 2 (x) - ^Pi(x) 4 , 

_ 30885 2 8669 4 1159 6 J^3_ 8 x 10 
~~256~ X + ^U X + l68 X + 2048 X + 512' 



1 

9~6 
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h{x) = P b {x) - Px{x)P A {x) - P 2 (x)P 3 (x) + P 1 (x)P 2 (x) 2 - P 1 (x) 3 P 2 (x) 

+P 1 (x) 2 P 3 (x) + ^P l (xf, 
5 

916731 2 33171 4 6453 6 823 8 319 10 7 12 

tJC 3> 3j 3j * 

512 128 256 512 5120 6144 ' 

h{x) = P 6 (x) - Pi(x)P 5 (x) - P 2 (x)P i (x) + Pi(x) 2 P 4 (x) - \p 3 {x) 2 + 2P 1 (x)P 2 (x)P 3 (x) 

-P 1 (x) 3 P 3 (x) + -U(*) 3 - Ip^xfP^xf + P 1 (xfP 2 {x) - \Pi(xf, 
6 2 6 

65518401 2 19425763 4 752825 fi 43783 8 3481 in 1255 12 

= x 2 H x 4 H x b H H x lu H x 12 

2048 4096 1536 4096 2048 24576 

and so on. ~ /3(^) coincide with the results in Ref. 10), where explicit 

expressions of f n (x) are given only for n < 3. It is interesting that the polynomials 
f n {x) are given in terms of P n {x) appearing in the naive perturbative expansion in a 
closed form. Tt should be instructive to remark here that /j(x)'s (i = 1,2,3, ...) are 
the cumulanwB of the sum 2~2^=o £ n Pn{x) in the sense that 

oo oo 

]T e n P n (x) ~ exp[^ e n Ux)\. (3.25) 

n=0 n=0 

In short, the RG method based on the construction of an envelope certainly 
resumes the perturbation series of the wave function and the resultant expression 
are given in terms of the cumulants of the naive perturbation series. Conversely, our 
method provides a mechanical way for computing cumulants of given series. 

3.3. Reproducing the WKB result 

We now examine how the WKB result Eq.(3.2) can be reproduced from the 
perturbation series obtained above. Bender and Bettencourt found that if all terms 
beyond l/512-e 4 a; 10 are neglected, the sum of the highest power terms in fj(x) (j < 4) 
is nicely rewritten as -rc 2 /4(l + 2ex 2 + 17/12 • e 2 x 4 + 5/12 • e 3 x 6 + 47/1152 • e 4 ^ 8 ) 1 / 8 , 
which behaves for large x as 

- v / e|a;| 3 /4(1152/47) 1/8 ~ y/e\x\ 3 /5.96G85, (3.26) 

in an excellent agreement with the WKB result. How about the higher orders. In the 
fifth order, the sum of the highest powers may be rewritten by neglecting all terms 
beyond 7e 5 x 10 /1286 as -x 2 /4-(l+ea; 2 /4-e 2 x 4 /24+e 3 x 6 /64-e 4 x 8 /128+7e 5 x 10 /1286) 

x 2 /4(l + 5/2 • ex 2 + 115/48 • e 2 x 4 + 35/32 • e 3 x 6 + 15/64 • e 4 x 8 + 4459/164608 • 

eV ) 1 / 10 . For large x, the coefficient of -^/e\x\ 3 is 4(164608/4459) 1 / 10 ~ 5.73827, 
which deviates from 6 more than the fourth order result. Unfortunately, the sixth 
order becomes worse: The sum of the highest powers is rewritten as — x 2 /4(l + 3ex 2 + 
29/8 • e 2 x 4 + 9/4 • e 3 x 6 + 577/768 • e 4 x 8 + 67621/493824 • e 5 x 10 + 1324349/35555328 • 
e 6 x 12 ) ly/12 , which makes the coefficient of — y^lxj 3 for large x 

4(35555328/1324349) 1/12 ~ 5.26181. (3.27) 
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This is actually plausible because the convergence radius of the perturbation series 
is zero; the cumulant series should be at best an asymptotic series. 

§4. Summary and Concluding Remarks 

In summary, we have successfully applied the RG methodic as formulated in 
Ref. |2|),[5]) to Schrodinger equation of the quantum anharmonic oscillator: The naive 
perturbation series of the wave function are resummed by the RG equation. We have 
found the following: Although the sum of the highest power in f n (x) can be or ga nized 
so that it becomes asymptotically proportional to y^M 3 as was done in , the 
coefficient of it reaches the closest value to 6, the WKB result, in the fourth order, 
then goes away monotonously from the closest value in the higher orders. We remark 
that the RG method as developed here can be also applied to the first excited state. 

It should be stressed that the method presented here can improve perturbative 
wave functions for all cases where naive perturbative solutions of Bender- Wu type are 
given; namely, the RG method combined with the Bender- Wu perturbation method 
constitutes a new powerful method for improving perturbative wave functions. It 
is interesting to extend the present method and apply it to quantum field theory, 
although it has been indicated that the notion of envelopes is useful for improving 
the effective potential given in the loop expansion, a* 

The present RG method is an application of the perturbative RG equation. 
Can nonperturbative RG equations be useful to construct global wave functions? 
It should be. Indeed a variational method called the delta-expansion methodic 
which also utilizes an RG- like equation but nonperturbatively has been extended 
for obtaining wave functions. The key ingredient of the extension is to construct 
an envelope of a set of perturbative wave functions as in the RG method but with 
a variational parameter. In this method, although the basic equation can not be 
solved analytically but only numerically, uniformly valid wave functions with correct 
asymptotic behavior are obtained in the first-order perturbation even for strong 
couplings and for excited states. In the present method, the basic equations are 
solved analytically, and the asymptotic form of the pSjave function is constructed 
explicitly, although a further resummation devised inED'B is needed for obtaining 
the asymptotic form. In this sense, the two methods are complementary to each 
other. It would be interesting if one can combine the two methods. 
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